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INTRODUCTION 


Bulk  polymeric  materials  are  especially  useful  in  aerospace  applica- 
tions because  of  their  high  strength  to  weight  ratio.  An  example  of  this  from 
the  fiber  industry  is  DuPont's  Kevlar.  In  order  to  "tailor" bulk  polymeric 
materials  for  specific  applications,  through  processing  or  initial  synthesis, 
knowledge  of  the  inherent  characteristics  of  the  material  is  required. 
Examining  each  individual  material's  morphology  and  correlating  this  with 
its  mechanical  properties  is  an  approach  to  ultimately  'tailor"  these  mater- 
ials. If  a reliable  (straightforward)  technique  existed  for  determining  the 
distribution  of  relaxation  times  associated  with  a given  mechanical  model, 
a direct  one-to-one  correlation  should  exist  between  this  distribution  and 
the  material's  morphology.  In  addition,  the  predictability  of  long  term 
mechanical  properties  would  be  available.  It  is  the  purpose  of  this  investi- 
gation to  present  in  an  illustrative  fashion  such  a technique  capable  of  yield- 
ing a reliable  distribution  of  mechanical  relaxation  times. 
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BACKGROUND 


The  problem  of  inferring  mechanical  relaxation  spectra  in  linear  vis- 
coelastic theory  is  herein  treated  as  a mathematically  ill-posed  problem. 

In  previous  investigations  the  established  regularization  technique  of  Tik- 
honov (Reference  1,2)  was  applied  to  the  mathematically  ill-posed  problem 
of  inferring  a reliable  molecular  weight  distribution  from  ultracentrifuge 
sedimentation  equilibrium  data  at  a single  angular  velocity  (References  3-5). 
This  technique  was  further  applied  (Reference  6)  to  computer  generated  data 
using  the  well  known  Fredholm  integrals  of  the  first  kind  encountered  in 
linear  viscoelastic  theory,  namely 


and 


E^(t)  = J H(T)  e(-t/r)  dlnT 


rw  2 2 
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(3) 


where  E (t)  is  a stress  relaxation  modulus,  E'(w)  and  E"(w)  are  the  in-phase 
r 

and  out-of-phase  components  of  the  complex  dynamic  modulus,  respectively, 
and  H(T)  is  the  sought  after  relaxation  spectra. 

The  difficulty  often  mentioned  in  finding  the  inverse  solution  to  equa- 
tions (1),  (2)  or  (3)  is  experimental  error  in  determining  E (t),  E'(w)  or 

r 

E"(w)  can  make  the  determination  of  H(T)  futile.  As  a result  iterative 
approximations  have  been  devised  so  as  to  control  and  limit  the  occurrence 
of  "wild  oscillations"  (References  6-12). 

The  present  technique  has  been  successfully  applied  to  computer  gen- 
erated analogous  experimental  data  (Reference  6).  Therefore,  it  was  felt 
that  application  to  real  experimental  data  was  now  necessary.  The  data 
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chosen  is  the  stress -relaxation  Master  Curve  data  of  NBS  polyisobutylene 
by  Tobolsky  and  Catsiff  (Reference  13).  In  the  following,  a brief  review  of 
the  theory  of  regularization  is  presented.  An  attempt  is  made  to  convey  this 
theory  in  light  of  the  method  of  undetermined  multipliers. 
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THEORY 


Rewrite  the  specific  equations  (1-3)  in  a more  general  form 

rh 

Q[x,f(y)]  = 0(x)  = / K(x,y)f(y)dy  (4) 

^ a 

where  0(x)  is  the  experimentally  determined  function,  f(y)  is  the  spectral 
response  function  and  K(x,y)  is  the  kernel  of  the  appropriate  integral  equa- 
tion. 

Assume  the  distribution  f(y)  is  continuous  and  bounded.  Therefore, 
let  a •&  y £ b and  similarly  c £ x £ d.  This  simply  establishes  upper  and 
lower  limits  on  the  observed  relaxation  times  and  frequencies  of  measure- 
ment, respectively.  Furthermore,  the  functions  0(x)  belong  to  some  com- 
plete metric  space  ^ and  the  functions  f(y)  to  another  complete  metric  space 
F,  i.e. , 0 € and  f e F.  Measuring  values  of  0(x)  experimentally  and  know- 
ing the  functional  form  of  the  integral  operator  which  formally  maps  f(y) 
into  0(x),  the  inverse  solution,  i.e.,  inferring  f(y)  given  j (>(x),  is  said  to  be 
a well  posed  problem  if:  (a)  the  solution  exists  for  f € F,  (b)  the  solution  is 
unique  in  F,  and  (c)  the  solution  depends  continuously  on  the  metrics  of  F 
and  ty.  Tikhonov  uses  a slightly  different  concept  of  correctness.  In  addi- 
tion to  the  spaces  and  F and  the  integral  operator,  he  assumes  a closed 
set  Ac  F,  f € A and  all  0 c A'.  A'  is  the  image  of  A after  application  to  the 
space  F of  the  integral  operator,  i.e.  , the  kernel  and  integration  process 
(Reference  14).  Instead  of  minimizing  the  conventional  functional 

N[f(y);0(x)]  -f  {Q[x,f(y)]  - 0(x)}2  dx  (5) 

J c 

Tikhonov  suggests  minimizing 

M“[f(y  );0(x)]  = N[f(y);0(x)]  + afl(n)[f(y)]  (6) 
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where  0 n is  the  regularizing  functional,  by  varying  tt  and  choosing  that  f(y) 

in  correspondence  with  the  inf  ||Ma  |[ . In  the  present  application 

n a=o 

/.b 


ar?^n)[f]  = E a.  f [f(l)(y)]2  dy 
i=o  ^ a 


(7) 


where  f^\y)  is  the  ith  derivative  of  f(y)  with  respect  to  y.  Now  each  O' is 

successively  varied  holding  all  others  constant.  In  the  present  application 

one  O',  was  sufficient, 
i 

From  the  method  of  undetermined  multipliers  presented  in  elementary 
calculus  (Reference  15),  a physical  picture  of  the  regularization  used  in  this 
investigation  can  be  achieved. 


Most  scientists  are  familiar  with  the  procedure  for  obtaining  the  max- 
ima and  minima  of  a function  on  the  assumption  that  the  variables  in  the 
function  are  linearly  independent.  However,  when  the  variables  are  con- 
nected by  certain  relationships,  i.  e.  , the  variables  are  not  all  linearly 
independent,  the  maxima  and  minima  are  conditional  extrema.  Consider  the 
function 


g(x  , X , . . .X  , X . , . . .X  ) 

1 2 m m+1  m+n 


of  (m+n)  variables  x.,  connected  by  n relationships 


K(xi,x2,...,xm,xm+i,...,xm+n)  = 0 (i=l>  2,  . . . , n) 


(8) 


(9) 


Let  the  function  g attain  a conditional  extrema  at  the  point  Q(x^,x^, 

. . . ,x  ).  Assuming  the  existence  of  the  derivatives  at  the  point  Q,  the 
total  differential  of  g must  vanish  at  Q, 

m+n 

S,  = 0 (10) 

s = l s 

At  the  same  point  Q,  the  following  n equations  hold  [differentiate  relation- 
ships  (9)1 
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m+n  5$ . 

£ r— — dx  = 0 (i=l,  2,  . . . , n) 
, ox  s 
S = 1 s 


(11) 


Multiplying  these  n boundary  functions  by  the  factors  X , \ ...  ,\  (Lagrange 

14  n 

multipliers),  and  adding  term  by  term  to  equation  (10),  one  obtains 


= 0 


(12) 


Usually  these  n factors  are  defined  such  that  the  n differentials  of  the  depend- 
ent variables  vanish. 


It  seems  reasonable  to  assume  that  the  present  integral  equations 
under  consideration  [equations  ( 1 ) -( 3 ) or  equation  (4)]  could  have  variables 
which  are  not  linearly  independent,  especially  when  "large"  errors  exist  in 
the  experimentally  determined  functions,  i.e.  {(b(x)).  These  errors  might 
not  seem  large,  but  if  the  inverse  transformation  needed  to  infer  jf(y)}  acts 
as  an  "amplifier"  of  noise.  Then  minimizing  N[f(y);0(x)]  i.e.,  eqn.  (5)  will 
not  yield  absolute  extrema  but  rather  conditional  extrema.  Not  knowing 
exactly  the  constraints  required  to  control  this  "noise"  one  would  be  led  to 
apply  the  logic  implied  by  equation  (12).  That  is,  assume  various  deriva- 
tives of  f(y)  are  the  constraints  and  vary  the  X's  (in  our  nomenclature  a.) 
until  a minimum  in  the  error  analysis  is  achieved.  The  "best  fit"  solution 
f*(y)  will  then  be  that  in  correspondence  with  the  smallest  or  inferior  of  all 
the  individual  conditional  minima.  It  is  just  this  technique  which  has  been 
applied  using  equations  (6)  and  (7)  along  with  the  computational  scheme  of 
quadratic  programming  (Reference  16).  The  unique  advantage  attained  in 
use  of  the  quadratic  programming  algorithm  is  that  only  positive  values  for 
f(y)  are  admissible.  Recall  f(y)  is  the  sought  after  spectral  response  func- 
tion and  negative  values  would  not  have  physical  meaning. 
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RESULTS 


The  present  technique  using  regularization  and  quadratic  programming 
was  applied  to  the  experimental  stress  relaxation.  Master  Curve,  data  pub- 
lished by  Tobolsky  and  Catsiff  (Reference  13),  The  values  of  log  E^(t)  are 
presented  in  Table  1 along  with  the  corresponding  logarithm  of  the  observa- 
tion times.  It  should  be  noted  that  even  though  the  data  appear  to  be  smooth 
on  a logarithmic  scale,  on  an  absolute  scale  values  may  have  large  fluctua- 
tions. That  is,  if 


[log  E (t.)] 

r x exp 


[log  E (t  )3  + 6 


(13) 


then 


where  6 = log  € 


E (t.)J  = e E (t.)L  = 106E  (t.)L  (14) 

r i 'exp  r i True  r i 'True  ' ’ 

Computationally  we  deal  with  actual  values  in  equation  (4).  These 
fluctuations  in  the  experimental  data  have  a profound  influence  on  the 
inverse  operational  matrix,  which  is  required  for  inferring  f(y).  The  results 
of  the  present  computation  are  numerically  given  in  Table  1.  The  relaxation 
spectrum  was  not  computed  at  the  same  relaxation  times  as  those  listed  in 
the  first  column.  Therefore,  the  inferred  log  H(T)  and  log  T values  are 
presented  in  the  last  two  columns.  They  are  spread  out  in  the  display  for 
ease  of  comparison  between  log  H(T)  and  log  H(T').  The  results  are  graph- 
ically displayed  in  Figures  1 and  2.  In  Figure  1 the  inferred  distribution  is 
represented  by  circles  and  the  'best  fit"  parameterized  "wedge  and  box"  dis- 
tribution of  Tobolsky  and  Catsiff  is  the  solid  line.  Using  the  circled  data 
of  Figure  1 or  the  data  in  the  last  two  columns  of  Table  1 and  equation  (3), 
the  computed  stress  relaxation  values,  column  four  of  Table  1 and  circled 
points  in  Figure  2,  were  determined.  To  achieve  these  results,  each 
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individual  f^\y)  [equation  (7)]  was  independently  applied.  The  correspond- 
ing OL  was  varied  and  the  distribution  in  correspondence  with  the  minimum 
of  coefficient  of  variation  between  the  experimental  E^  ^gft)  anc^  ca^” 

culated  E (t)  was  used.  The  chosen  distribution  was  in  correspondence  with 
r 

the  inferior  of  these  minima  coefficients  of  variation.  In  the  present  case 
for  i=l. 


was  the  regularization  function  yielding  inf||E  (t)  -E  (t)  ||. 

r exp  r calc 


(15) 
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CONCLUSION 


A computational  technique  is  of  value  to  an  experimenter  if  it  is  reli- 
able, easy  to  use,  and  does  not  consume  a large  amount  of  high  speed  digital 
computer  time.  In  the  present  case,  the  reproducibility  of  the  log  H(T)  dis- 
tribution almost  identical  to  that  achieved  by  the  na  priori11  functional  'Wedge- 
box’1  shape  is  very  encouraging.  The  only  restriction  on  the  sought  after 
distribution,  in  the  present  case,  is  that  the  boundaries  of  the  range  of 
relaxation  times  be  broad  enough  so  that  the  function  and  its  derivatives 
vanish  at  the  boundaries  and  that  the  function  is  continuous.  The  present 
technique  is  easy  to  use  because  of  the  simplicity  of  the  mathematical  con- 
cepts. Finally,  the  computational  time  required  for  a broad  search  of 
numerical  values  e.  g.  , 10  <:  oc  ^ 10^  required  only  about  one  second 

central  processing  time  on  a CDC  CYBER  70  computer.  With  experience 
an  investigator  should  know  which  to  search  and  depending  upon  his  usual 
experimental  error,  the  approximate  range  of  numerical  values  to  scan. 

This  will  further  reduce  the  actual  computational  time. 

Having  the  capability  to  easily  infer  a relaxation  spectra  from  experi- 
mental mechanical  data  within  the  regime  of  linear  viscoelasticity,  one  is 
led  to  apply  the  technique  to  data  from  a variety  of  specimens.  As  mentioned 
within  the  text  this  will  allow  a correlation  between  mechanical  model  relax- 
ation spectra  H (T)  and  the  observed  morphology.  This  should  be  a step 
toward  eventual  'tailoring11  molecular  level  processes  to  yield  desired 
macroscopic  end  use  morphologies.  It  is  realized  that  this  is  a long  term 
process.  However,  another  more  immediate  use  of  this  new  capability  will 
be  the  time  saving  in  being  able  to  do  a few  experiments  and,  using  existing 
viscoelastic  theories,  predict  other  unmeasured  mechanical  behavior  of  the 
material.  This  will  especially  be  of  value  for  long  time  mechanical  behavior. 
Again  one  should  remember  that  the  predictions  are  only  as  good  as  the 
theories  used  in  the  extrapolation  of  information.  That  is,  there  will  always 
be  the  need  for  experimental  checking,  but  the  number  of  experiments  per- 
formed will  decrease  as  confidence  in  the  theory  increases. 
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Finally,  since  the  present  investigation  concerns  itself  only  with  linear 
viscoelastic  theory,  the  greatest  impact  will  come  by  extending  these  ideas 
to  the  nonlinear  viscoelastic  regime. 
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Figure  2.  The  Stress  Relaxation  Curves  in  Correspondence  with 
the  Two  Relaxation  Spectra  Presented  in  Figure  1. 


TABLE  1 


Comparison  of  Mechanical  Relaxation  Spectrum  When  "a  priori" 
Wedge -Box  Shape  is  Assumed3,  and  When  No  Previous 
Functional  Form  is  Assumed  (Regularization  Technique  Applied) 

* 2 

(The  stress -relaxation  data  (dyne/cm  ) calculated  is  also  presented) 
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(log  T' ) 

hr 

10.  52 
10.  50 
10.49 
10.47 
10.  46 
10.45 
10.43 

9.  81 

-14.2 

10.41 
10.  38 
10.35 

9.  82 

-13.0 

10.  31 

9.  87 

-12.  4 

10.27 
10.21 
10.  15 
10.07 
10.01 

10.  10 

-11.  8 

9.  89 

10.11 

-11.2 

9.79 

9.68 

9.  55 

9.  81 

-10.7 

9.40 

9.24 

9.41 

-10.  1 

9.11 

i 

8.99 

8.  87 

8.98 

-9.5 

8.76 

8.  59 

8.45 

8.  69 

-8.9 

8.32 
8.  19 
8.01 
7.94 
7.  81 

8.  15 

-8.  3 
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TABLE  1 (Continued) 
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TABLE  1 (Concluded) 


log  t 
log  T 
(hr) 

Tobolsky  & Catsiff3- 

Regularization  Technique 

log 

Er,298(t) 

log  H 
(log  T) 

log 

Er(t)calc 

log  H 
(log  T') 

log  T' 
hr 

-0.  8 

6.  59 

6.22 

6.  59 

.6 

6.  55 

6.22 

6.  55 

6.  32 

-0.6 

.4 

6.  50 

6.22 

6.  52 

.2 

6.45 

6.22 

6.  46 

.0 

6.39 

6.22 

6.40 

6.24 

0.0 

0.2 

6.33 

6.22 

6.  34 

.4 

6.26 

6.22 

6.26 

.6 

6. 17 

6.22 

6. 18 

6.27 

0.6 

.8 

6.08 

6. 19 

6.09 

1.0 

5.97 

6.14 

6.00 

.2 

5.  85 

6.08 

5.  87 

6.  13 

1.2 

.4 

5.70 

6.00 

5.72 

.6 

5.  54 

5.  89 

5.  53 

. 8 

5.34 

5.74 

5.  33 

5.  82 

1.  8 

2.0 

5.  18 

5.  57 

5.  14 

.2 

4.90 

5.38 

4.94 

.4 

4.  50 

5.  18 

4.71 

5.  18 

2.4 

.6 

4.00 

4.98 

See  reference  13. 
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